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ABSTRACT 

It has been suggested that chondrules and calcium-aluminum-rich inclusions 
(CAIs) were formed at the inner edge of the protoplanetary disk and then en- 
trained in magnetocentrifugal X-winds. We study trajectories of such solid bodies 
with the consideration of the central star gravity, the protoplanetary disk grav- 
ity, and the gas drag of the wind. The efficiency of the gas drag depends on a 
parameter n, which is the product of the solid body size and density. We find 
that the gravity of the protoplanetary disk has a non-negligible effect on the tra- 
jectories. If a solid body re-enters the flared disk, the re-entering radius depends 
on the stellar magnetic dipole moment, the disk's gravity, the parameter rj, and 
the initial launching angle. The disk's gravity can make the re-entering radius 
lower by up to 30%. We find a threshold r], denoted as T) t , for any particular 
configuration of the X-wind, below which the solid bodies will be expelled from 
the planetary system. t] t sensitively depends on the initial launching angle, and 
also depends on the mass of the disk. Only the solid bodies with a rj larger than 
but very close to rjt can be launched to a re-entering radius larger than 1 AU. 
This size-sorting effect may explain why chondrules come with a narrow range of 
sizes within each chondritic class. In general, the size distributions of CAIs and 
chondrules in chondrites can be determined from the initial size distribution as 
well as the distribution over the initial launching angle. 

Subject headings: celestial mechanics — stars: protostars — protoplanetary disks 
- meteorites, meteors, meteoroids — comets: general — minor planets, aster- 
oids: general 

1. Introduction 

X-winds are high-speed bipolar collimated jets around young stellar objects (YSOs) 
powered by enhanced magnetic activities and disk-magnetosphere interactions (Shu et al. 
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1996, 1997). An X-wind can be driven magnetocentrifugally from the inner edge of the disk 
where accreting gas is either diverted onto stellar field lines to flow onto the star or to be 
flung outwards with the wind (Shu et al. 1994a,b). 

X-winds are also suggested to be the underlying mechanism for the production of 
calcium-aluminum-rich inclusions (CAIs) and chondrules for our Solar system (Shu et al. 
2001). The subsolar reconnection ring (X-region) at the inner edge of the protoplanetary 
disk provides high-temperature environment to form CAIs and chondrules, and the X-winds 
can lift the processed rocks back to colder regions of the disk. At least 50% by mass of all 
the magnesium-iron silicates that finally constructed planetesimals and planets may have 
been processed as chondrule or CAIs (Tayler 1992; Shu et al. 1996). 

Isotopic measurements also provide evidence of the X-wind origin of planetary materi- 
als. CAIs have 18 0/ 16 and 17 0/ 16 ratios that are several percent lower than any other 
meteoritic and planetary materials in the Solar system (Clayton 1993). This isotopic feature 
can be explained by photochemical self-shielding of carbon monoxide in the solar nebula, 
which requires that planetary materials passed through the inner part of the accreting pro- 
toplanetary disk (Clayton 2002). Again the processed materials need to be returned back to 
colder regions by X-winds. 

However, there are also other formation mechanisms that have been proposed to explain 
the chemical, mineral and isotopical properties of chondrules, such as exothermic chemical 
reactions, nebular lightning, magnetic reconnection flares, gas dynamic shock waves in the 
protoplanetary nebula, radiative heating, etc. (see Rubin 2000 for an extensive review). 
Desch et al. (2010) critically examines the X-wind model for the formation of chondrules 
and CAIs against isotopical and mineralogical evidences, but also advocates detail modeling 
on the fate of the particles launched by the outflow to see if the particles can be retained in 
proper annuli of the disk. 

The parent bodies of chondritic meteorites originate mostly from the asteroid belt at 
about 2.5 astronomical units (AU) from the Sun. Recent observations find evidence of 
CAI-type and chondrulelike particles present in a short-term comet 81P/Wild 2 (Swindle & 
Campins 2004; Brownlee et al. 2006; Flynn et al. 2006; Nakamura et al. 2008). Wooden 
(2008) concluded from mineralogy investigations that these refractory grains should be pre- 
accretionary with respect to the formation of asteroids and efficient outward radial transport 
mechanism is required. This discovery requires large-scale radial transport mechanisms, such 
as X-winds and the turbulent diffusion of particles (e.g. Ciesla 2009). However, it has not 
been studied whether X-winds are able to drive solid bodies to the radius of the Kuiper belt, 
or >30 AU. 
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It has not yet been clearly explained under what conditions X-wind can deliver chon- 
drules and CAIs to the asteroid belt or the Kuiper belt. Shu et al. (1996) and Shang (1998) 
provides several trajectories of solid bodies in X-winds in the dimensionless form without 
the consideration of protoplanetary disk gravity nor the vertical geometry of the disk. Liff- 
man (2005) investigated the size-sorting effect of chondrules from magnetic pressure driven 
jet flows, but only for distance <3AU from the Sun. In this paper we study the transport 
of solid bodies (e.g. CAIs and chondrules) by the X-winds, with the protoplanetary disk's 
gravity and vertical geometry taken into account. We also investigate whether and where 
the launched solid bodies re-enter the disk and how the re-entering radius depends on the 
stellar magnetic dipole moment and the size and density of the solid bodies. 



2. Model 

2.1. X-wind Configurations 

In this section we provide an introduction on the X-wind configurations based on Shu et 
al. (1994a, 1996), and identify several stellar and disk parameters that govern the X-wind. 

An accreting protoplanetary disk that interacts with the rotating magnetosphere of a 
protostar with a pure dipole is truncated at a radius of 

where /x* is the stellar dipole moment, M* is the stellar mass, Md is the inflow rate of the 
disk, G is the gravitational constant, and C\ is a dimensionless coefficient that we assume 
to be unity in the following (e.g. Shu et al. 1996). In steady state, the inner disk edge 
corotates with the star at the Keplerian angular velocity of 

The wind has a terminal velocity v w when leaving the X region. The terminal velocity 
can be written as 

v w = (2J g - 3) 1 / 3 i? x Q x , (3) 

where J g is the dimensionless specific angular momentum of the gas (Shu et al. 1994). 

The formation of CAIs requires the thermal treatment at 1800~2200 K. CAIs and 
chondrules entrained in X-wind experience the peak temperature when they are lifted into 
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the direct sunlight. The peak temperature is given by Shu et al. (1996) as 

f r /R*\ /R*\/ Rl\ 1/2 l > \ 1/i 

^^iT^+s^i^H^J-UJl 1 -^; J) ' (4) 

where is the stellar luminosity, er is the Stefan-Boltzmann constant, and i?* is the radius 
of the star. 

From above we see that the physical properties of X-winds are determined by three 
parameters of a planetary system: M*, and M-q. In theory, the X-wind can exist in either 
the 'embedded' stage, when the protostar is still embedded in its natal envelope of gas and 
dust; or the 'revealed' stage, when the outflowing wind has reversed the infall of the envelope 
and revealed the central star and the accretion disk (Shu et al. 1996). At different stages, 
the star has different masses, the stellar magnetic dipole moment varies in different ranges, 
and the accretion rate has different values. As long as the lifetime of these two stages are 
fairly well constrained by the observations of young stellar objects, we use typical values of 
M* and Md for each stages, and let /x* vary in different ranges, respectively. Following Shu 
et al. (1997), we define a high state, an average state and a low state according the value of 
/i* for each stage, as shown in Tabled) Key X-wind properties Rx, fix, v w and T peak of each 
state are also given in Table [TJ As evident from Table (U a protosun with stronger magnetic 
field truncates the protoplanetary disk at a larger radius, and therefore has lower T pea k and 
v w . Notably, CAIs can only form during the EA and EL cases. 

For each state, the gas density and flow velocity of the X-wind can be determined 
numerically, from a set of equations involving the conservation of mass, specific momentum 
and specific energy (see Shang 1998 for the detail formulation). In the following, we formulate 
the equation of motion in the inertial frame, with the cylindrical coordinates (I, 0, z) centered 
at the star and with z = defining the midplane of the disk. The asymptotic behavior of 
the X-wind dynamics sufficiently far from the X-region (i.e. I > 50Rx) is: (1) the poloidal 
velocity reaches the terminal velocity; (2) the toroidal velocity decreases as / _1 ; (3) the 
density of the wind falls off as l~ 2 . Assuming the streamlines to be straight lines originating 
from the X-region and the open angle of the bipolar wind to be 8 W = 45°, we follow the 
procedure of Shang (1998) to construct the X-wind as the background for the motion of the 
solid body. 

2.2. Extended Model 

In this section we set up the framework to calculate the trajectories of solid bodies 
entrained in the X-wind. A solid body with size of R s and density of p s entrained in a 
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X-wind is subject to three forces: the gravity of the star, the gravity of the protoplanetary 
disk, and the aerodynamic drag of the X-wind. The total force exerted on the solid body is 

,, GM*m s 

-r - A(l,z)m s e z - B(l,z)m s ei + F diag , (5) 



where m s is the mass of the solid body, and r is the distance from the solid body to the center 
of star. The second and the third terms in the right-hand side of equation (jSJ) represent the 
protoplanetary disk gravity, which has not been included in previous trajectory calculations. 
A(l, z) and B(l, z) can be calculated by the superposition of gravitational contributions from 
surface elements throughout the disk, viz. 

= _d_ /*- r 2£g£2fW£ (6) 

dz J Rin J v ///2 + /2 + 2 2 _2/'/cos0' ' 

dl J Rin Jo ^/l' 2 + l 2 + z2-2l'lcos(j)> ' 

where I' and <$' are dummy variables of integration spanning the whole disk, R m and i? ou t 
are the inner and outer edge of the protoplanetary disk, and E(Z) is the surface density of 
the disk. In the calculation of disk's gravity, we assume the geometry of the disk to be a 
thin plane at z — 0, and use R lVL = Rx and R out = 40 AU. The integrals can be evaluate 
numerically once we know the surface density profile of the protoplanetary disk. 

We adopt that the disk has the mass distribution as the standard minimum mass solar 
nebula (MMSN) given by Hayashi (1981). However, the density of the protoplanetary disk 
might be much higher. For example, Desch (2007) estimated a denser MMSN based on the 
NICE model. The disk's surface density in our model is parameterized as 

£(/)=/xl700(— ) gem" 2 , (8) 

where / is a multiplying factor specifying the mass of the disk and the rest is the MMSN 
profile given by Hayashi (1981). Note that Desch (2007) suggested the surface density 
exponent to be —2.2 instead of —1.5, but in this work we adopt the profile by Hayashi 
(1981). The acceleration due to the disk gravity A(l,z) and B(l,z) can then be computed 
numerically. 

The protoplanetary disk has the vertical extension that may intercept the solid body 
trajectory and influence the radius where the solid body re-enters the disk. Assuming the 
disk is supported by the thermal pressure, the scale height of the disk is 
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where k& is the Bolzmann constant, p is the mean molecular mass in the atomic mass unit 
(for the solar abundance p = 2.3), M H is the mass of hydrogen atom, and T(l) is assumed 
to be the temperature of the MMSN disk as T(l) = 280 l~ 1/2 K (see Hartmann 2009). In 
the following we regard the scale height in equation (jH]) as the 'boundary' of the disk. Once 
the vertical position of the solid body becomes smaller than the scale height of the disk, 
z < H, the solid body re-enters the disk, and then the trajectory calculation is terminated. 
We do not intend to follow the trajectories of the solid bodies inside the disk. Once the 
horizontal position of the solid body becomes outside the disk, or I > 40 AU, the solid body 
is ejected from the planetary system. In this work, we assume a thin disk when calculating 
the gravitational force from the disk. Still, we do consider the vertical geometry of the disk 
using a realistic "flared disk" model when determining the radius where the solids re-enter 
the disk. We will demonstrate that the vertical structure of the disk has a significant effect 
on the re-entering radius. 

The Epstein gas drag applies when the relative velocity of the solid body and the wind 
is less than the sound speed of the wind, and the radius of the solid body is smaller than 
the mean free path of the gas. For example, the mean free path of the wind is 10 ~ 10 2 cm 
during the embedded stage, and 10 2 ~ 10 3 cm during the revealed stage. Therefore when 
calculating the trajectories of chondrule-or-CAI-size solid bodies entrained in X-winds, the 
Epstein gas drag law should be used, viz. 

F drag = ~7ri&U(v«, - v s ) 2 + c 2 ] 1 / 2 ^ - v s ) , (10) 

where p w , c w , and v w are the density, the sound speed and the bulk velocity of the wind, 
and v s is the velocity of the solid body (Weidenschilling 1977; Shang 1998; Youdin & Chiang 
2004). The gas drag depends on the density and velocity of the wind. The acceleration due 
to the gas drag is 

a drag = ^ = — [K - v s ) 2 + c 2 ] 1 / 2 (v w - v.) , (11) 

Tflg Tj 

where we define the drag efficiency parameter rj = R s p s . The gas drag is more effective for 
a solid body with smaller r], meaning smaller or less dense. 

Assuming that a solid body is launched at the initial angle of #0 with respect to the 
midplane of the disk from the reconnection ring at the speed of the local X-wind, we integrate 
the equations of motion (jSJ) numerically to find the trajectory of the solid body. The allowing 
range of 9 is [0, 9 W }. In summary, the model calls 6 parameters: 3 parameter to specify 
the X-wind, M„, Md and p*; / to specify the mass of the protoplanetary disk; 77 to specify 
the efficiency of the gas drag; and #0 to specify the initial launching angle of the solid body. 
For any particular state in Table [U the former 3 wind-related parameters are fixed, and the 
latter 3 parameters are free to be explored. 
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3. Results 

3.1. Trajectories of Solid Bodies 

The trajectories of solid bodies of different size and density (or 77) in the EA and RA 
states are plotted in Figure [TJ Here we assume / = 5 and only explore the situations of 
maximum lift (9q = 9 W ) and intermediate lift (6q — 9 w /2). It is evident from Figure [1] that 
smaller and less dense objects can be launched to a larger orbit, in agreement with previous 
works (e.g. Shu et al. 1996). Note that although the left and the right panels of Figured] 
appear to be very similar, they show trajectories of solid bodies with very different 77 (see 
the caption for details). Based on Figured] the X-wind in the EA state can lift much larger 
and denser (or larger 77) solid bodies to a significantly large radius than that in the RA state 
does. 

The vertical extent of the disk intercepts some trajectories. From Figure [I] the trajec- 
tories approach straight lines in the I — z plot at large I and z, in agreement with previous 
works of Shu et al. (1996), Shang (1998) and Liffman (2005). Without the consideration of 
disk's flared vertical geometry, these straight-line trajectories may easily be interpreted as 
'ejection'. For example, in the EA case, / = 5 and 9 = 9 W , the solid body with rj = 1.3 
cgs unit re-enters the disk at I = 8.5 AU (see the left panel of Figure [1]). If the disk was 
assumed to be a thin surface, such object would be ejected from the system. Therefore, one 
has to consider the disk's vertical geometry in order to yield realistic diagnostics on whether 
the solid bodies re-enter the disk or at what radius the solid bodies re-enter the disk. 

The initial launching angle #0 has a large effect on the trajectories of the solid bodies. 
From Figure [H a solid body is launched to a much larger orbit when 9o — 9 W . For example, 
in the EA case and / = 5, a solid body with 77 = 1 cgs unit is ejected from the system if 
#0 — 9 W , but re-enters the disk at ~ 2 AU if 9 Q = 9 w /2 (see Figured]). For another example, 
in the EA case and / = 5, a solid body with rj = 1.3 cgs unit re-enters the disk at ~ 8.5 AU 
if 6*o = W , but re-enters the disk at much smaller radius (~ 0.7 AU) if #0 = 9 w /2. Therefore, 
if numbers of solid bodies of such 77 are initially launched from various angles in [0, 9 W ], the 
X-wind spreads them to the disk's various annuli up to ~8.5 AU. In general, under the same 
condition, a solid body can be lifted to a larger orbit and thus re-enters the disk at a larger 
radius if initially launched from a larger angle. 

There exists a threshold value of rj of solid bodies for any particular X-wind configura- 
tion, disk's mass and initial launching angle, denoted as r] t , below which the solid bodies will 
be expelled from the planetary system (see Figured]). This result is expected because the 
aerodynamic drag on the light particles by the wind is strong, i.e., these particles are more 
strongly coupled to the gas. For the heavy (77 > r] t ) particles, on the other hand, drag forces 
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are much less important compared to the gravitational forces from the central star and the 
disk; and as a result, they fall back into the disk. 

3.2. The Effect of Disk Gravity 

In this work we include the gravity of the protoplanetary disk into the calculation of 
trajectories of solid bodies entrained in X-winds for the first time. Here we explore the effect 
of the disk gravity, and demonstrate that the disk mass modifies the trajectories and makes 
the re-entering radius lower significantly. 

In Figure [2] we plot the trajectories of the same solid body with different disk mass 
factor /. If / = 0, the disk's gravity is not considered and we return to the situation of 
Shu et al. (1996). From Figure [2j we find that the disk's gravity has a negligible effect 
on the trajectory at / < 1 AU. Detail analysis on the forces reveals that at I < 1 AU, 
the gravity of the central star provides the major force in the — e z direction. When I > 1 
AU, the trajectories under different / become different. At such distance, the gravity of 
the protoplanetary disk becomes comparable with the gravity of the central star in the —e z 
direction. As a result, the solid body flies closer to the disk (smaller z at the same I) when 
the disk is more massive (larger /). Therefore, the re-entering radius becomes also smaller 
when the disk is more massive. Compared the MMSN disk case to the non-disk case, the 
re-entering radius decreases by ~ 10%. If still considered the fact that actual protoplanetary 
disk is very likely to be much more massive than the MMSN disk, for example / = 5, the 
re-entering radius can be lowered by up to ~ 30%. 

Here we find that the protoplanetary disk provides a non-negligible force in the — e z 
direction and modifies the trajectories at I > 1 AU. The previous calculations of Shu et 
al. (1996) that did not consider the disk gravity are therefore imperfect in estimating the 
re-entering radius. 

3.3. The Retention of Solids 

In this section we study the sensitivity of r] t on the disk's mass / and the initial launching 
angle 9q. We show how rjt depends on / and 9q i n Figure [3j By definition, solid bodies of 
rj < rjt entrained in X-winds are expelled from the Solar system. As mentioned in Section 
2.2, we consider a solid body to have been expelled if / > 40 AU in the trajectory calculation. 
Note that the value of r] t will change if this cutoff changes. 

From Figure [3j the threshold rjt depends most sensitively on the initial launching angle 
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6q. Typically when / < 10, rjt at the maximum launching angle 9q = 9 W is one-order-of- 
magnitude greater than that at a low launching angle of 9$ = 0.39 w . For example, during 
the EA state and / = 5, % ~ 0.98 cgs unit at 9 = 9 W ; and r) t ~ 0.09 cgs unit at 9 = 0.39 w . 
Similarly, during the RA state and / = 5, r) t ~ 0.039 cgs unit at 9 = 6 W \ and r\ t ~ 0.006 cgs 
unit at 9 = 0.39 w . 

When the disk becomes more massive, the solids need to be smaller and lighter in order 
to be expelled by the X-winds. As seen in Figure El r] t decreases by ~10% when / increases 
from to 10, and rj t decreases by ~30% when / increases from to 100. The dependence of 
rjt on / is less sensitive, but should also be taken into account when considering the retention 
of solids in the protoplanetary disk. 

The X-wind is stronger during the embedded stage, as a result, rj t during the embedded 
stage is much larger than that during the revealed stage, as shown in Figure[3J It is interesting 
to note that T) t during the revealed stage is compatible with the typical chondrule size range 
in meteorites, if the density is assumed to be 2 ~ 4 g cm -3 . Finally, an important observation 
that we can make from Figure [T] and Figure [3] is that the horizontal distance that a solid 
body can travel depends on its size and density (77) very sensitively. 



3.4. The Size-Sorting Effect 

In this section we investigate how the re-entering radius depends on the stellar magnetic 
dipole moment and size and density of the solid bodies. We calculate the trajectories of solid 
bodies of different rj in the states listed in Table [TJ and plot the relations between 77 and 
the re-entering radius for each state in Figure HI There are still flexibility in choosing the 
value of parameter / and 9q. In a particular state, a solid body can be launched to the 
maximum re-entering radius when / = and 9q = 9 W , as plotted in Figure HI There is no 
theoretical minimum of the re-entering radius for any particular state, since 9$ can be as low 
as 0. However, we still plot the re-entering radius with / = 5 and #0 = 9 w /2 on Figure IU to 
indicate how small the re-entering radius can be for the intermediate lift of any particular 
state. Again, it is evident that smaller bodies re-enter the disk at larger radii. 

The main result that one can draw from Figure S] is that only the solid bodies whose 
77 are larger than but very close to r\ t of a particular state can be launched to the annuli of 
asteroid belt or Kuiper belt. For example, during the EH state and / = 0, 9 — 9 W , r\ t ~ 0.7 
cgs unit, and the solid bodies with 77 ~ 1.1 - 1.3 cgs unit re-enter the disk at the asteroid-belt 
radius (see the solid curve of the top panel of Figure H]). For another example, during the 
RH state and / = 0, 9 = 9 W , rj t ~ 0.03 cgs unit, and the solid bodies with 77 ~ 0.04 - 
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0.05 cgs unit re-enter the disk at the asteroid-belt radius (see the solid curve of the bottom 
panel of Figure H]). The densities of chondrules and CAIs are in a narrow range around 2-4 
g cm -3 (e.g. Hudges 1980). As a result, rj t is in general equivalent to the threshold size. Of 
course, meteoritic parent bodies may have a wider feeding zone than the current asteroid 
belt location during their formation. However, the result does not change even if we consider 
a wider feeding zone, because in order to be lifted to a radius larger than 1 AU, the solid 
body must already have a size close to the threshold size (see Figure HJ). 

In any particular state, the location of the asteroid belt, where we envisage the parent 
bodies of carbonaceous chondrites, can only receive solid bodies with sizes larger than but 
very close to the threshold size. However, the variation of initial launching angle 6q affects r\t 
significantly (see Figure [3]) and enlarges the acceptable range of r\ significantly. In reality, we 
know very little about the distribution of 9 . The X-wind model (Shu et al. 2001) suggests 
that the magneto-rotational instability (MRI) happens at the X-point to transport the rocks 
from the reconnection ring to the launching height. Based on the turbulent nature of the 
process, it is reasonable to assume that the distribution of 8$ is very wide, or nearly uniform. 
Imagine the CAIs and chondrules form in the X-region and are launched in X-wind from a 
wide range of 9 . Due to this variation of 9 , in any particular state, the annuli of the asteroid 
belt can receive a wide spectrum of solid bodies with different rj, in other words, different 
sizes. For example, during the EH state, if the condition (/ = 0, 6> = 6 W ) is specified, the 
range for rj to re-enter the asteroid-belt radius is 1.1~1.3 cgs unit (see the solid curves of the 
top panel of Figure H]). However, once the variation of #0 from 9 W to 9 w /2 is considered, the 
range for 77 to re-enter the asteroid-belt radius becomes as wide as 0.6~1.3 cgs unit. In the 
previous section we find that r] t depends most sensitively on 6q, so here the enlargement of 
the t] range is mainly due to the the variation in 60 ■ 

Liffman (2005) predicts very effective size sorting and very sharp size distributions in 
a constant outflow. This is because in Liffman (2005), the initial angle of articles is fixed, 
corresponding to a single curve in our Figure HI Here we step forward to suggest that the 
variation of 6q may make the size distribution of CAIs or chondrules in chondrites wider. 

The second result is that an annulus of the protoplanetary disk receives one-order- 
of-magnitude smaller solid bodies (assuming the density is generally the same) during the 
revealed stage than during the embedded stage. For example the location of asteroid belt 
receives solid bodies of 77 ~ 0.8-1.7 cgs unit during the EA state, but solid bodies of 77 ~ 
0.04-0.07 cgs unit during the RA state, with the uncertainty of / and 9 (~ [6 W /2,9 W ]) 
considered. Table [TJshows that only the EA and EL states can reach the required temperature 
Tpeak > 1800 K and produce CAIs. In contrary, chondrules need lower temperature and are 
younger than CAIs (e.g. Scott 2007). It is then likely that chondrules mainly formed during 
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the revealed stage. If this is the case, the chondrules should be smaller than CAIs in any 
particular chondritic group because of the X-wind transport. However, such difference could 
not be found in chondrites (e.g. Scott 2007). We will address to this discrepancy in the next 
section. 

Last but not least, the variations of re-entering radius according to //* are different 
for small and large bodies. As the star becomes less magnetized, the small solid body re- 
enters the disk at a larger radius, but the large solid body re-enters the disk at a smaller 
radius. This feature may counter the physical intuition. Actually, as //* decreases, or the 
star becomes less magnetized, the accretion disk is truncated at a smaller radius. However, 
the star rotates much faster according to equation ([2]) so that the terminal velocity of X- 
wind becomes faster. For small bodies, the wind terminal velocity determines how much 
momentum it can gain during the aerodynamic launching phase; therefore, less magnetized 
star produces faster winds, and then solid bodies re-enter the disk at larger radii. For large 
bodies, their aerodynamic launching phase is very short and their re-entering radii are close 
to the truncation radius Rx- 



4. Discussion 

In any particular state of X-winds, /, #0 and rj determine whether an entrained solid 
body re-enters the disk, as well as the re-entering radius. Section 3.3 shows that the retention 
of solid bodies depends most sensitively on the initial launching angle 9q. Especially, small 
bodies fall back to the disk if they are launched at a sufficiently small launching angle. In 
order to estimate the percentage of entrained solid bodies that eventually fall back, one has 
to know the distribution of 6q, which relates to complicated processes of retrieving materials 
from the reconnection ring. 

Although blurred by the uncertainty in 8q, the size sorting effect of X-winds may have 
significant implications in our understanding of CAI and chondrule populations. Because of 
the uncertainty in 8q, the present model can only provide an upper limit of rj below which a 
solid body can be delivered by X-winds to certain location of the disk. 

Based on Figure HJ the acceptable 77 for a solid body to be delivered to the location of 
the asteroid belt during the revealed stage is up to 0.08 cgs unit. If we assume chondrules 
mainly formed during the revealed stage and the density of chondrules is in the range of 
2^4 g cm -3 (e.g. Hudges 1980), the X-wind transport model here predicts the maximum 
chondrule radius to be ~ 0.04 cm, compatible with all carbonaceous chondrite groups except 
the CV class (e.g. Wurm & Krauss, 2006). In average, CV chondrites have largest chondrules 
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with size of ~ 0.1 cm. From Figure HI chondrules in CV chondrites cannot be delivered to 
the location of the asteroid belt during the revealed stage. Instead, they might form during 
the EH state, or during the transitional stage between the embedded stage and the revealed 
stage. 

In particular, CH chondrites have smallest chondrules with size of ~ 0.002 cm (Jones 
et al. 2000). According to our X-wind transport model, a solid body as small as the 
chondrules in CH chondrites can easily be expelled from the system. One possible way to 
retain such tiny objects is to have very small initial launching angle 9$. However, if the 
chondrules that formed in the X- region were launched isotropically at all angles in [0, 9 W ], 
large ones should be delivered simultaneously to the location of the asteroid belt. Therefore 
chondrules need to be launched preferably at small launching angles in order to explain the 
depletion of large chondrules in CH chondrites. Another possible explanation is that the 
chondrules of CH chondrites were delivered at the very end of the revealed stage, when the 
disk accretion and the X-wind approached to the end. If this is the case, chondrules in CH 
chondrite experienced the thermal treatment of a lower temperature, which could be tested 
by mineralogical examinations. 

CAI-type and chondrulelike particles found in a short-term comet 81P/Wild 2 are very 
small having size of ~ 20 /im (e.g. Brownlee et al. 2006; Nakamura et al. 2008). X- winds 
could easily expel these tiny objects from the Solar system (see Figure Hj). However, the 
particles had been larger than the present size because they were ablated and disaggregated 
during the deceleration in the aerogel (Flynn et al. 2006). The entrance- hole sizes of impact 
tracks place an upper limit for the size of the initial incoming particles, which is about 0.1 
cm (or r] ~ 0.3 cgs unit). From Figure IU the solid bodies with such sizes can be transported 
to the Kuiper belt during the embedded stages. 

CAI sizes are also different in various chondritic groups (Hezel et al. 2008), correlating 
with the chondrule sizes. CV chondrites appear to have largest CAIs with size range in 
0.01~0.15 cm, and CH chondrites have smallest CAIs with size generally <0.01 cm (e.g. 
Zhang & Hsu 2009). Because CAIs formation requires high temperature, only the EA and 
EL states are suitable (see dotted and dashed curve on the top panel of Figure 0J. We find 
that CAIs in all chondritic groups are lighter than the threshold r] t of 9q — 9 w /2 during the 
EA and EL states. Therefore, for these tiny primitive grains to be retained in the disk, the 
current X-wind model requires that they are launched preferably at low initial launching 
angles (9 < 9 w /2). This requirement provides a new aspect to test the X-wind model, and 
also suggests the importance to understand the initial 9 distribution. 

Finally, according to Figure H] large solid bodies whose r\ well exceed rj t cannot be lifted 
over long distances. They will re-enter the inner disk and follow the accreting inflow. It 
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is possible that they enter the reconnection ring for the second time, during the timescale 
of embedded or revealed stage. After being thermally processed for the second time, they 
are expelled again by X-winds and then may enter a larger orbit if they experience some 
fragmentation. Especially, independent enveloping compound chondrules have a secondary 
layer of chondrulelike material enclosing an unrelated primary, indicating that they experi- 
ence multiple melting (Wasson et al. 1995). It would be very interesting to study the role 
of X-wind in the multiple thermal treatments of such chondrules. 

5. Conclusion 

In this work we calculate the trajectories of solid body entrained in X-winds with the 
consideration of the protoplanetary disk's gravity and geometry. We find that the disk gravity 
becomes comparable to the gravity of the central star at I > 1 AU, and the disk's mass affects 
the trajectories and the location where the solid bodies re-enter the disk significantly. The 
gravity of a disk whose mass is 5 times of the MMSN disk can lower the re-entering radius 
by up to 30%. 

The product of the size and the density, denoted as i] in this paper, is the key parameter 
to evaluate the efficiency the gas drag and the sorting effect. We find a threshold value of 
i], denoted as r) t , for any particular X-wind state, disk's mass and initial launching angle. 
Solid bodies with rj < rj t are expelled by the X-wind from the Solar system. We find that 
rjt depends on the initial launching angle 6q most sensitively, but also depends on the disk's 
mass factor /. 

We demonstrate again the size-sorting effect of the X-wind transport by showing that 
only the solid bodies with a r\ larger than but very close to r) t can be launched to the location 
of the asteroid belt. We also find that the variation of the initial launching angle may smooth 
out the size-sorting effect, because solid bodies with a wider range of 77 can be delivered to the 
desired location. This improved X-wind transport model is consistent with the chondrule 
size in all carbonaceous chondrite groups other than CV and CH. The size of chondrules 
in CV and CH chondrites, as well as the size of CAIs in carbonaceous chondrites, require 
additional mechanisms to be explained. 

In order to compare the model prediction with the measured chondrule size distribution, 
it is necessary to further calculate the i] distribution after the X-wind size-sorting. To do 
so, one has to know the initial size distribution as well as the distribution over the initial 
launching angle, which are determined by the formation mechanism of chondrules in the 
X-region. Furthermore, there is little constrain on the variation of the stellar magnetic 



- 14 - 



dipole moment. The final r\ distribution reflects the entire evolution history of the early Sun. 
For example, if the stellar activity remains stable for a long period of time, we should see 
relatively sharp size distributions of chondrules. Therefore, the 77 distribution of chondrules 
may contain valuable information of the early Solar System. 

R. Hu is grateful to S. Seager for helpful advice in the preparation of the manuscript, to 
L. Elkins-Tanton for the encouragement of publication of this work, to J. Meyer and B. Weiss 
for helpful discussions, and to the anonymous referee for the improvement of the manuscript. 
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Table 1: Definition and X-wind properties of all states considered in this paper. Each state is 
defined by a set of parameters (M*, Mb, fi*)- For each case, key X-wind properties Rx, &x, 
v w and Tpeak are calculated based on equations (jHHJ). Full names of cases are given following 
the table, but simplified names, given in the first column of the table, are used in the text. 
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-L peak 
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2 For EH, EA and EL cases, L* = 4.4L Q and R* = 3R & are used in the calculation of 
Tpeak, where L & and R Q is the current luminosity and radius of the Sun. Similarly, for 
RH, RA and RL cases, we assume = 2.5L and R* = 3R & . 
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Fig. 1. — Trajectories of solid bodies entrained in X- winds of EA state (left) and RA state 
(right) with / = 5. In both panels the solid lines are the trajectories of maximum initial 
launching angle (9q = 9 W ), and the dashed lines are the trajectories of intermediate initial 
launching angle (9 = 9 w /2). The dotted lines show the scale height of the disk. Both panels 
have the same vertical axis, so the tags are only shown in the left panel. In the down-right 
region of each panel is the enlarged figure showing the detail of trajectories in I < 2 AU. 
Left: From top to bottom the solid lines (or dashed lines) correspond to trajectories of rj ~ 
0.6, 1.0, 1.3 and 1.6 cgs unit, respectively. Note that for 9 = 9 W , the solid bodies of r\ ~ 
0.6 and 1.0 cgs unit are expelled from the planetary system; and for 9q = 9 w /2, none of 
these solid bodies is expelled from the planetary system. Right: From top to bottom the 
solid lines (or dashed lines) correspond to trajectories of rj ~ 0.01, 0.03, 0.05 and 0.06 cgs 
unit, respectively. Note that for 6 = 9 W , the solid bodies of r\ ~ 0.01 and 0.03 cgs unit are 
expelled from the planetary system; and for 9 = 9 w /2, only the solid body of r\ ~ 0.01 cgs 
unit is expelled from the planetary system. 
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Embedded Stage - Average State 



ri=1.0 cgs unit 




l[AU] 

Fig. 2. — Top: Trajectories of a solid body of 77 = 1.0 cgs unit entrained in the X-wind in the 
EA state with #0 = 9 w /2. Bottom: Trajectories of a solid body of r\ = 0.04 cgs unit entrained 
in the X-wind in the RA state with 8 = 6 w /2. In both panels the trajectories correspond 
to the disk mass factor / = 0, 1, 3, 5, respectively. The dotted lines show the scale height 
of the disk. Note that the vertical and horizontal scales in these panels are different. 
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Fig. 3. — The sensitivity of the threshold r] t on / and # in the EA state (top) and the RA 
state (bottom). For each line, we fix the value of # as shown on the figure, vary / and find 
the corresponding rj t . 
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Fig. 4. — Relations between the solid body size and density (rj) and the re-entering radius. 
Vertical light dashed lines indicates the radii of the asteroid belt (AB) and the Kuiper belt 
(KB) of the Solar system. For both stages, the solid curves describe the high state, the 
dashed curves describe the average state, and the dotted curves describe the low state. Each 
state is described with two curves. The higher curve corresponds to / = and 9 = 9 W , and 
the lower curve corresponds to / = 5 and 9 = d w /2. 



